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ABSTRACT 

The  occurrence  of  the  El  Nino  Southern  Oscillation  (ENSO)  has  been  studied 
from  numerous  air-ocean  interaction  aspects,  yielding  theories  implying  a 
positive-only  correlation  between  the  time  rate  of  change  of  sea  surface 
temperature  anomaly  (dTs/dt)  and  the  corresponding  change  in  the  depth  of  the 
ocean  mixed  layer  (h').  However,  an  alternate  proposal  by  Chu  (1991a,b,  1992) 
holds  that  there  are  both  positive  and  negative  correlations  between  dTs/dt  and 
h'.  Chu's  proposal,  based  upon  the  modified  Kraus-Turner  ocean  mixed  layer 
(OML)  thermodynamic  model,  goes  further  to  say  that  the  regions  occupied  by 
these  positive  and  negative  correlation  values  correspond  roughly  to  those 
occupied  by  positive  and  negative  values  of  a  surface  forcing  function,  P,  that  is 
an  indication  of  the  strength  of  the  counteracting  surface  forcings  of  wind  work 
and  buoyant  damping.  The  model-generated  fields  of  surface  wind  stress,  net 
heat  flux,  sea  surface  temperature,  and  primary  (mixed)  layer  depth  are  analyzed 
for  a  two  year  ENSO  event,  1986  -  87,  covering  the  region  of  the  equatorial 
Pacific,  ION  to  10S  and  to  HOE  to  80W,  in  order  to  determine  the  validity  of 
these  theories.  The  results  shows  a  fairly  uniform  P-field  over  the  entire  period, 
consisting  of  a  negative  P-field  (indicative  of  a  shallowing  regime)  in  most  of  the 
western  Pacific  and  in  the  far  eastern  part  of  the  eastern  Pacific  with  a  positive  P- 
field  (indicative  of  an  entrainment  regime)  in  between.  The  correlations  between 
dTs/dt  and  h'  show  that  there  exists  both  positive  and  negative  correlations 
between  dTs/dt  and  h',  however,  there  is  no  appreciable  similarity  between  the 
regions  occupied  by  positive  and  negative  correlations  and  those  occupied  by 
positive  and  negative  P-values. 
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I.     INTRODUCTION 

A.    RATIONALE  FOR  STUDY 

Currently,  there  are  several  significant  theories  containing  simple  ocean 
thermodynamics  that  have  been  developed  to  study  the  phenomenon  of  the  El 
Nino  Southern  Oscillation  (ENSO).  Certain  theories  differ,  however,  concerning 
the  correlation  between  the  time  rate  of  change  of  sea  surface  temperature 
anomaly  (3T's/3t)  and  the  anomaly  of  upper  layer  thickness  (T|),  or  mixed  layer 
depth  anomaly  (h').  The  relationship  between  T|  and  h'  will  be  explained  in  the 
following  background  section,  but  for  now,  it  is  sufficient  to  say  that  they  are 
interchangeable. 

Numerous  published  theories  (Wyrtki,  1975;  McCreary,  1983;  Cane  and 
Zebiak,  1985;  Battisti,  1988,1989;  and  Hirst,  1986  among  others)  show  a  positive- 
only  correlation  between  dTjdt  and  T|.  The  modified  Kraus-Turner  (K-T)  ocean 
mixed  layer  (OML)  thermodynamic  model  proposed  by  Chu  (1991a,b,  1992) 
holds,  conversely,  that  there  are  both  positive  and  negative  correlations  between 
these  two  parameters  depending  on  the  physical  process  being  considered.  The 
purpose  of  this  study  is  to  analyze  data  to  establish  the  proper  correlation 
between  these  parameters  and  thus  establish  a  correct  basis  for  the  study  of  the  El 
Nino  phenomenon. 

L      Upwelling  Physics 

The  theories  based  on  upwelling  physics  are  repesented  by  the  follwing 
thermodynamic  equation  (Hirst,  1986): 

dTs/dt  +  u*Tx=k(ari-Ts),k>Oy(T>0  (1) 


which  implies  a  positive  correlation  between  3Ty3t  and  r|'.  In  the  limit  of  large  k, 
this  equation  reduces  to  T's  =  ar|,  indicating  a  positive-only  correlation  between 
TsandT|. 

2.      Ocean  Mixed  Layer  Physics 

Chu  (1991a,b,  1992)  draws  a  distinction  between  "upwelling  physics" 
and  "OML  physics"  in  dealing  with  the  correlation  between  dTjdt  and  h'.  When 
considering  "upwelling  physics,"  based  upon  a  cooling  upper  layer  where  the 
outgoing  heat  flux  from  the  base  of  the  mixed  layer  is  greater  than  the  incoming 
heat  flux,  the  following  relationship  is  valid: 

drsldt~h\  (2) 

where  there  is  a  positive  correlation  between  3T's  /  di  and  h\  Based  upon  this 
type  of  physics,  equation  (1)  of  the  Hirst  model  is  valid.  However,  when 
considering  "OML  physics",  where  the  upper  layer  may  be  under  either  strong  or 
weak  wind  forcing,  there  exist  both  positive  and  negative  correlations  between 
these  two  parameters 

dTsldt~{+l-)h\  (3) 

where  there  exists  a  positive  correlation  for  strong  wind  forcing  and  negative 
correlation  for  weak  wind  forcing.  Chu  takes  this  relationship  one  step  further, 
relating  these  positive  and  negative  correlations  to  corresponding  positive  and 
negative  values  of  a  surface  forcing  function,  P.  This  function  will  be  discussed 
further  in  the  following  section. 


B.  BACKGROUND 

1.      Relationship  of  Mixed  Layer  Depth  to  Upper  Layer  Thickness 
Anomalies 

The  value  of  the  upper  layer  thickness  anomaly,  r),  is  based  on  the 
position  of  the  main  thermocline,  and  is  defined  as: 

V  =  P/g'-d  (4) 

where  g'  =  gAp  /  p,  r|  is  the  instantaneous  change  in  the  depth  of  the  thermocline, 
and  d  is  the  initial  depth  of  the  mixed  layer.  Figure  1  illustrates  the  physical 
relationship  between  h  and  r\. 


Figure  1.  Illustration  of  upper  layer  anomaly,  r\,  where  d  represents 
the  initial  thickness,  and  h  represents  the  instantaneous  mixed  layer 
depth  (from  Yamagata,  1985). 


Essentially,  r\  represents  the  change  in  the  depth  of  the  thermocline 
lying  just  below  the  mixed  layer,  and  therefore  corresponds  to  the  mixed  layer 
depth  anomaly,  h'. 

2.      Significance  of  Surface  Forcing  function,  P 

A  discussion  of  P  should  begin  with  the  modified  K-T  model  that 
introduced  P  as  a  surface  forcing  function.    The  value  of  P,  whether  positive  or 


negative,  is  indicative  of  the  predominant  surface  forcing,  wind  work  or  buoyant 
damping  respectively,  on  the  OML,  and  is  defined  by  the  following  equation: 

P  =  Cxu\  -  C2aghwQQ  I  pw0cpw  (5) 

(a)  (b) 

where  Ci  and  C2  are  tuning  coefficients,  uw*  is  the  water  surface  friction  velocity 
defined  as  (CDPao/pwo)1/2  I  Vb  I ,  Vb  is  the  surface  wind  velocity,  hw  is  the 
thickness  of  the  OML,  Q0  is  the  net  surface  heat  flux,  paU  and  pwo  are 
characteristic  values  of  surface  air  density  and  water  density,  and  a  is  the  thermal 
expansion  coefficient.  The  values  used  for  the  tuning  coefficients,  Ci  and  C2,  are 
1 .0  and  0.5  respectively  (Schneider  and  Muller,  1 990).  Term  (a)  represents  the 
surface  forcing  due  to  wind  work,  and  term  (b)  defines  the  surface  forcing  due  to 
buoyant  forcing/damping.  Wind  work  acts  to  increase  turbulence  in  the  mixed 
layer,  and  thus  increase  the  depth  of  the  mixed  layer.  Buoyant  damping  is  caused 
by  downward  heat  flux,  +Q0,  which  acts  to  stabilize  the  mixed  layer  by  heating 
the  upper  portion  of  the  mixed  layer  leading  to  stratification  and  thus  a 
shallowing  of  the  mixed  layer.  Strong  surface  wind  forcing,  uw*,  will  result  in  a 
positive  value  for  P.  Likewise,  strong  downward  heat  flux,  Q0,  accompanied  by 
moderate  to  weak  wind  forcing,  will  result  in  a  negative  value  for  P.  Generally, 
the  value  of  P  varies  from  positive  to  negative  throughout  the  ocean  basin  based 
on  the  diurnal  and  seasonal  variations  in  QQ  and  the  synoptic  and  seasonal 
variations  in  uw*. 

The  modified  K-T  model  incorporates  this  P-value  in  the  following 
equation: 

ar  ,  idt + aw  •  vr,  =  Qo-K(P)Q-h  (6) 


'w 


PwOCpw"M> 

(a)  (b)  (c) 


where  term  (a)  is  the  rate  of  change  for  Ts,  term  (b)  the  horizontal  advection  of  Ts, 
and  term  (c)  the  thermodynamic  forcing,  including  the  heat  flux  at  the  base  of  the 
mixed  layer,  Q„h.  Of  significance  in  term  (c)  is  the  heavyside  step  function,  AW(P). 
This  function  is  referred  to  as  an  OML  "switcher."  For  positive  values  of  P,  the 
value  of  the  switcher  is  1,  and  for  zero  or  negative  values  of  P,  the  value  of  the 
switcher  is  0.  These  two  cases  have  the  following  effect  on  the  modified  K-T 
model  in  equation  (6): 
a.    P<0 

In  this  case,  the  value  of  A(P)  is  0,  yielding  the  following  equation: 
dTs/dt  +  uw.VTs= ^-—  (7) 

PwOcpwflw 

Here,  a  negative  P-value  results  from  weak  wind  forcing  and/or  strong  downward 
heat  flux,  causing  a  shallowing  in  the  OML.  Since  the  mixed  layer  will  heat  more 
rapidly  as  its  depth  decreases,  there  will  be  a  positive  correlation  between  h'  and 
3Ts3t.  This  case  is  illustrated  in  figure  2. 


Qo 


Weak  Surface  Wind  Forcing 
z  =  0 

Ocean  Mixed  Layer 
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Figure  2.  Kraus-Turner  OML  Thermodynamics  under  weak  wind 
forcing,  P<0  (from  Chu,  1991a). 

b.    P>0 

In  this  case,  the  value  of  A(P)  is  1,  yielding: 
dT f  Idt  +  uw  •  VT t  =  &Z%  (8) 


'  wQCpw'hv 


As  noted  before,  a  positive  P-value  results  from  strong  wind  forcing,  which  causes 
entrainment  of  deep  water  into  the  OML,  as  indicated  by  the  presence  of  the  Q.h 
term.  Since  Q.h  is  greater  than  QQ,  the  OML  is  referred  to  as  a  "cooling  layer" 
(Miller,  1976),  and  an  increase  in  the  mixed  layer  depth  will  be  accompanied  by  a 
corresponding  increase  in  dTs/dt.  Thus,  there  would  be  a  positive  correlation 
between  h'  and  9Ty9t.  In  physical  terms,  the  thicker  the  layer  becomes  the  less 
rapidly  it  will  cool.  This  case  is  illustrated  in  figure  3. 


Qo 
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Figure  3.  Kraus-Turner  OML  thermodynamics  under  strong  wind 
forcing,  P>0  (from  Chu,  1991a). 


II.  DATA  COLLECTION  AND  PROCESSING 

A.    DATA  COLLECTION 

In  order  to  effectively  study  the  correlation  between  dTs/dt  and  h'  and  its 
relationship  to  the  P- value,  the  following  data  sets  are  required:  monthly  mean 
net  surface  heat  flux  (QnetX  monthly  mean  surface  wind  stress  (tw),  daily  mean  sea 
surface  temperature(Ts),  and  daily  mean  mixed  layer  depth  (h).  The  area  of 
coverage  includes  the  majority  of  the  equatorial  Pacific.  Therefore,  the  study  area 
is  chosen  as  140  E  to  80  W  and  ION  to  10S  inclusive.  Figure  4  is  a  map  of  the 
equatorial  Pacific  region.  The  time  period  is  chosen  to  span  a  designated  ENSO 
event,  January  1986  -  December  1987.  Data  was  collected  for  this  study  from 
two  sources,  Fleet  Numerical  Oceanography  Center  (FNOC)  in  Monterey, 
California,  and  the  National  Center  for  Atmospheric  Research  (NCAR)  in  Boulder, 
Colorado.  The  data  products  received  from  FNOC,  Ts  and  h,  were  generated  by 
the  Thermodynamic  Ocean  Prediction  System  (TOPS)  -  coupled  Expanded  Ocean 
Thermal  Structure  (EOTS)  Analysis  model.  The  data  products  received  from 
NCAR  were  Qnet,  Tw,  and  Ts,  and  were  archived  from  the  European  Center  for 
Medium  Range  Weather  Forecasting  (ECMWF)  synoptic  forecasts.  A 
description  of  these  data  sets  and  their  forecasting  and  modeling  methods  is 
provided  below. 

L      FNOC  Archive  Data  Sets 

Data  was  acquired  through  FNOC  from  the  archives  of  the  Data  Base 
Management  Division.  During  the  time  period  of  1986-87,  the  data  products  of 
sea  surface  temperature  and  primary  (mixed)  layer  depth  were  generated  by  the 
TOPS-  coupled  EOTS  (TEOTS)  analysis  technique.   A  discussion  of  the  FNOC 


data  will  include  a  section  on  TEOTS  analysis  methods  and  another  section  on 
TEOTS  data  product  output. 
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Figure  4.  Equatorial  Pacific  region 
a.     TEO TS  A  nalysis  Meth od 

The  EOTS  analysis  scheme  combines  ocean  thermal  climatology 
with  real  time  thermal  observations,  applying  a  series  of  corrections  via  the  Fields- 
by-Information-Blending  (FIB)  analysis  technique.  The  ocean  climatology  is 
drawn  from  an  archive  of  coarse  resolution  (320  km)  monthly  climatology.  The 
observations  are  provided  by  real  time  ship  and  satellite  data.  Daily  Ship 
observations  consist  of  about  150-200  XBT's  ,  and  1200-2000  ship  injection/ 
bucket  temperature  measurements.  Daily  Satellite  measurements  consist  of 
50,000-80,000  Multi-Channel  Sea  Surface  Temperature  (MCSST)  retrievals  from 
NOAA.  This  system  is  limited,  however,  since  local  bursts  of  atmospheric  forcing, 
which  cause  considerable  change  in  the  mixed  layer  depth  and  sea  surface 
temperature  are  often  neglected  by  using  such  a  sparse  number  of  real  time 
measurements.  For  this  reason,  FNOC  adapted  a  system,  endorsed  by  Elsberry 
and  Garwood  (1980)  and  Clancy  and  Martin  (1981),  coupling  a  synoptic  mixed 
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layer  model  (TOPS)  to  the  basic  EOTS  analysis  to  generate  a  continuously 
updated,  real-time,  analysis-forecast-analysis  system.  The  gist  of  the  system  is  to 
initialize  the  TOPS  model  daily  using  the  EOTS  analysis  fields,  produce  a  forecast 
field  driven  by  observed  and  forecasted  atmospheric  forcing,  and  subsequently 
feed  this  field  back  into  the  EOTS  analysis  as  a  first-guess  for  the  following 
analysis  day.  This  system  is  collectively  referred  to  as  TEOTS,  and  appears  in  the 
form  of  a  flow  chart  in  figure  5  (FNOC,  1987). 
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OCEAN  THERMAL 

ANALYSIS/PREDICTION 

(TEOTS/TOPS) 

The  TEOTS  analysis  (nowcast)  provides 
initial  conditions  for  the  TOPS  prediction 
(forecast).  Thus,  TEOTS  updates  TOPS 
with  the  real-time  data.  The  TOPS 
prediction,  which  accounts  for  the 
response  of  the  surface  mixed  layer  to 
atmospheric  forcing,  is  combined  with 
climatology  and  fed  back  into  the 
subsequent  TEOTS  analysis  as  the 
first-guess  field  for  the  analysis. 


THERMAL  INPUT  TO  FNOC  ACOUSTIC 
MODELS 

GRAPHICS  PRODUCTS  TO  COMMAND  AND 
CONTROL  CENTERS 


Figure  5.  TEOTS  Flow  Chart  (from  FNOC,  1987) 

One  concern  in  deciding  whether  or  not  to  use  the  TEOTS- 

generated  data  involves  the  relationship  between  Ts  and  h.  One  of  main  goals  of 

this  study  is  to  determine  a  correlation  between  these  two  parameters,  and 

therefore  it  is  necessary  to  establish  that  there  are  no  physical  assumptions 


concerning  this  correlation  within  the  TOPS  forecast  that  might  bias  the  results  of 
this  study.  The  following  is  a  brief  discussion  of  the  physics  used  in  TOPS. 

The  TOPS  consists  of  four  prognostic  equations  representing  the 
conservation  of  temperature,  salinity,  and  momentum  in  the  upper  ocean  (FNOC, 
1987). 


dt      dz  dz       pwc  dz      ax 

d  ,     ~.      d         —  d  T     d  T 


dS       d  .    — r-=r        dS.      d  .     7;.      d  ,     -=,. 

—  =  — (-w  S'  +  v — ) (uaS)-  —  (vaS) 

dt      dz  dz       dx  dy 

d(     Tx     Atd2S     d2S^ 


(9) 


(10) 


du        -      d       — r-7        du.         - 

-r-  =  /v+  —  (-w'u'  +  v—)-Du  (11) 

at  oz  oz 


"3-,  ~\  ~\ 

?l  =  -fZ  +  —  (-^r7+v— )-Dv  (12) 

dt  dz  dz 


where  T  is  the  temperature,  S  is  the  salinity,  u,  v,  and  w  are  the  x,  y,  and  z 
components  of  the  current  velocity,  F  is  the  solar  radiation  flux  (positive 
downward),  pw  is  the  reference  density,  c  is  the  specific  heat  of  seawater,  D  is  a 
damping  coefficient,  v  is  a  diffusion  coefficient,  z  is  the  vertical  coordinate 
(positive  upward  from  the  surface),  and  ua,  va,  and  wa  are  the  x,  y,  and  z 
components  of  the  "advection"  current.    The  parameter  A  is  a  horizontal  eddy 
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diffusion  coefficient  whose  value  depends  on  the  horizontal  grid  spacing. 
Equations  (9)  -  (12)  exclude  the  effects  of  geostrophic  flow,  which  was 
determined  to  have  a  negligible  effect  on  mixed  layer  dynamics. 

These  equations  were  parameterized  by  the  following  boundary 
conditions  (Mellor  and  Durbin,  1975  and  Clancy  and  Martin,  1981): 

df 


WT=-tqSH 


dz 


3  r 

w>S'=-£qSH—- 
dz 


du 


Wu'=-lqSM  — 


dv 


WV=-ZqSM  — 
dz 


(13) 
(14) 
(15) 
(16) 


where  £  is  the  turbulence  length  scale,  q  is  the  square  of  twice  the  turbulent 
kinetic  energy  and  is  defined  by  the  following  equation: 


WM 


ydz) 


+ 


\dz) 


+  tqSH 


(  g_dp 


Py,      dZ 


\ 


15^ 


=  0 


(17) 


The  three  terms  in  equation  (17)  represent  wind  forcing,  or  shear  production, 
buoyant  damping,  and  viscous  dissipation  respectively.  The  variables  Sm  and  Sh 


are  functions  of  the  gradient  Richardson  number,  Ri,  defined  as  follows: 

8   dp 


Ri  = 


Pw  & 


(du\2    (dv\2 


ydz  j 


Adz 


(18) 
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The  Richardson  number  basically  indicates  the  predominant  forcing  between 

shear  production  and  buoyant  damping.    If  the  Richardson  number  exceeds  a 

critical  number,  Ri  =  .23,  buoyant  damping  stratifies  the  mixed  layer  sufficiently 

enough  to  inhibit  the  turbulent  mixing  of  shear  production,  and  the  layer  enters  a 

shallowing  regime  (Miles,  1961  and  LeBlond  and  Mysak,  1978).  Conversely,  if  Ri 

<  .23,  shear  production  dominates,  resulting  in  entrainment  and  a  deepening  of 

the  mixed  layer.   The  extent  to  which  the  mixed  layer  deepens  is  dependent  on 
the  magnitude  of  the  turbulent  length  scale,  £y  defined  as  follows: 

I     \z\qdz 
£  =  0A±z (19) 

As  seen  in  the  progression  through  the  prognostic  equations  (9)  -  (12)  to 
the  closure  of  equation  (19),  there  are  no  assumptions  made  regarding  the  direct 
relationship  of  mixed  layer  depth  to  sea  surface  temperature.  The  TOPS  model  is  a 
very  meticulous  mixed  layer  forecasting  system  whose  only  assumption  excludes 
the  presence  of  geostrophic  flow.  For  this  reason,  the  TEOTS  analysis  system  was 
determined  to  be  unbiased  regarding  the  objective  of  this  study. 

b.     TEOTS  Analysis  Output  Description 

The  output  of  the  TEOTS  analysis  took  the  form  of  two  different 
grid  schemes,  one  for  Ts  and  one  for  h.  Mixed  layer  depth  is  output  on  a  grid 
using  a  polar  stereographic  projection.  This  type  of  scheme  consists  of  63  X  63 
grids  overlaid  on  maps  of  the  Northern  and  Southern  Hemispheres  whose  origins 
are  the  poles.  The  grid  spacing  is  approximately  320  kilometers.  A  disadvantage 
of  this  grid  scheme,  when  dealing  with  equatorial  data,  is  that  the  resolution  of  the 
data  decreases  radially  outward  from  the  polar  origin,  yielding  the  worst 
resolution  at  the  equator.    Sea  surface  temperature  is  output  on  a  grid  using  a 
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spherical,  or  Mercator,  plot.  This  type  of  scheme  consists  of  a  73  X  144  grid 
overlaid  on  a  Mercator  projection  of  the  globe.  The  grid  spacing  for  the  spherical 
plot  is  approximately  2.5°,  or  280  kilometers.  The  contoured  plots  of  monthly 
mean  mixed  layer  depth  are  found  in  Appendix  D.  The  contoured  plots  of 
monthly  sea  surface  temperature  are  found  in  Appendix  C. 
2.      NCAR/ECMWF  Archived  Data  Sets 

As  previously  mentioned,  the  ECMWF  data  sets  were  received  from 
NCAR.  A  description  of  the  numerical  models  used  to  generate  this  data  will  be 
quite  brief  due  to  the  emphasis  on  mixed  layer  dynamics  rather  than  the  physics 
of  the  processes  driving  atmospheric  forcing. 

a.     ECMWF  Numerical  Weather  Prediction 

The  wind  and  surface  flux  fields  provided  by  ECMWF  are 
generated  by  an  operational  model  using  numerical  weather  prediction.  The 
model  used  for  wind  analysis  and  prediction  produces  6-hourly  forecasts  on  15 
standard  levels  in  the  vertical  and  on  a  1.125°  horizontal  grid.  The  model 
combines  observations  obtained  within  a  six  hour  window  of  the  analysis  time 
with  a  multivariate  height,  wind,  and  thickness  system  (Lorenc,  1981)  from  the 
preceding  analysis  period  to  produce  a  forecast  for  the  following  period.  The 
method  used  to  combine  these  two  inputs,  statistical  interpolation,  assigns 
weights  to  the  these  inputs,  and  uses  geostrophy  and  other  error  checking  criteria 
to  smooth  the  results  (Lonnberg,  1982).  A  system  flow  chart  is  shown  in  figure  8. 
Once  the  surface  wind  field  has  been  generated,  the  surface  wind 
stress  fields  are  computed  using  the  relation: 

r0=pCDU2  (20) 
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where  p  is  the  density  of  air  in  the  planetary  boundary  layer,  Cd  is  the  drag 
coefficient  over  the  surface,  and  U  is  the  10  meter  wind  provided  by  the 
foregoing  model. 
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CURRENT  OBSERVATIONS 

-Surface  land  and  sea  reports 
-Radiosonde  and  pilot  reports 
-Satellite  thickness  reports 
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STATISTICAL  INTERPOLATION 

Combines  previous  forecast  and 
current  observations,  applying 
the  linear  constraint  of 
geostrophy  in  order  to  provide 
consistent  analyses.  Also 
includes  error  detecting  and 
correcting  schemes. 
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|NEW  FORECAST  FIELDS 

Figure  6.  ECMWF  analysis  /  forecast  flowchart 

Surface  flux  fields  are  simply  a  byproduct  of  the  standard  outputs  (winds, 
temperature,  potential  temperature,  and  humidity)  of  the  ECMWF  model. 
Shortwave  and  longwave  radiation  are  computed  using  straightforward 
atmospheric  physics  combined  with  satellite  cloud  imagery  to  form  an  algorithm 
modeling  the  behavior  of  the  radiation  budget.  Radiation  affects  the  ocean 
mixed  layer  both  directly,  by  increasing  or  decreasing  the  heat  flux  driving  the 
buoyant  damping  on  the  layer,  and  indirectly,  by  affecting  the  development  of 
clouds  and  weather  patterns  which  have  a  local  effect  on  shear  production.  A 
schematic  of  these  radiative  processes  is  shown  in  figure  9.  Longwave  radiation 
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emitted  by  the  surface  is  computed  using  the  modified  Stefan-Boltzman  law 


RL=-€oT 


(21) 


where  e  is  the  emissivity  of  water,  a  is  the  Stefan-Boltzman  constant,  and  T  is  the 
temperature  of  the  surface  in  degrees  K.  Since  longwave  radiation  is  emitted  not 
only  by  the  surface  but  by  clouds  and  the  atmosphere,  the  estimation  of 
longwave  radiation  is  a  fairly  complicated  function  of  the  earth  and  satellite- 
based  soundings.  Shortwave  radiation  at  the  top  of  the  atmosphere  is  computed 
using  an  inverse  square  law  coupled  with  Lambert's  cosine  law  to  yield 

RQ=S0(dm/d)2cosy  (22) 

where  So  is  the  solar  constant,  dm/d  is  the  ratio  of  the  mean  to  the 
actual  distance  from  the  sun,  and  y  is  the  solar  zenith  angle.  In  order  to  obtain  the 
actual  surface  value  of  shortwave  radiation,  absorption  and  scattering 
predominantly  by  clouds  are  also  considered. 


n'.tr.uv, 


short    waves 


long    waves 


Figure  7.  Radiative  Processes  in  the  Planetary  Boundary  Layer  (from 
Wergen,  1982) 


Latent  heat  flux  is  computed  from  the  basic  equation 
HL  =  LeE 


(23) 
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where  Le  is  the  latent  heat  of  evaporation/condensation  and  E  is  the  rate  of 

evaporation/condensation  determined  by  the  soundings.    Sensible  heat  flux  is 

calculated  using  the  bulk  aerodynamic  approach  based  on  the  Monin-Obukhov 

similarity  theory.    This  theory  assumes  that  near  the  surface,  vertical  wind  and 

potential  temperature  gradients  depend  only  upon  height,  surface  sensible  heat 

flux,  surface  momentum  flux,  and  the  thermal  expansion  coefficient  of  air. 

Accordingly,  ECMWF  uses  the  temperature  and  winds  at  their  lowest  sounding 

level  (33  m)  and  the  surface  temperature  to  compute  the  sensible  heat  flux  at  the 

surface  as 

^  =  [i  /i*#     J2*G(h/z0,Ri)\Vh\A6h  (24) 

ln(/z/z0) 

where  k  is  von  Karman's  constant,  h/zo  is  the  ratio  of  height  to  roughness  length, 

Vh  is  the  wind  velocity  at  the  33  m  level,  A0h  is  the  change  in  potential 

temperature  between  the  surface  and  the  33  m  level,  and  G(h/zo,Ri)  is  an 

analytical  function  that  is  continuously  updated  based  on  observations  and 

performance  of  the  model  (Wergen,  1982).  The  Richardson  number,  Ri,  is  given 

by 

Ri  =  ^t  (25) 


6\Vh\ 

and  is  an  indication  of  stability  within  the  surface  boundary  layer. 
b.     ECMWF  Output  Description 

The  ECMWF  data  fields  accessed  through  NCAR  were  archived  in 
grid-in-binary  (GRIB)  format.  In  order  to  compute  the  surface  friction  velocity, 
the  following  data  sets  were  extracted:  u-component  of  surface  wind  stress,  and 
v-component  of  surface  wind  stress.  In  order  to  compute  the  net  surface  heat 
flux,  the  following  data  sets  were  extracted:  net  surface  shortwave  radiation,  net 
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surface  longwave  radiation,  surface  sensible  heat  flux,  and  surface  latent  heat 
flux.  The  data  fields  are  derived  from  6  hour  forecasts  valid  for  0  Z  and  12  Z,  and 
are  plotted  on  a  1.125°  grid.  The  contoured  plots  of  monthly  mean  net  surface 
heat  flux  are  found  in  Appendix  A.  The  contoured  plots  of  monthly  mean  water 
surface  friction  velocity  are  found  in  Appendix  B. 

B.    DATA  PROCESSING 

Data  processing  consisted  of  two  basic  steps.  First,  wind  and  heat  flux  data 
fields  are  combined  with  mixed  layer  data  to  compute  values  for  the  surface 
forcing  function,  P.  Sea  surface  temperature  and  mixed  layer  depth  were  then 
analyzed  to  determine  the  correlation  between  3Ts/3t  and  h'. 

1.      P- Value  Computation 

The  P-value  fields  were  computed  by  equation  (4).  Since  the  majority  of 
the  data  required  for  this  computation  came  from  the  ECMWF  forecasts,  a  1.125° 
grid  was  chosen  to  represent  the  P-field.  One  P-field  was  calculated  per  month 
from  January  1986  -  December  1987.  The  following  equation  was  used  to 
compute  the  value  of  the  net  surface  heat  flux,  Qo: 

Qo  =  Qsw  +  Qtw  +  Qs  +  Qi  (26) 

where  Qsw  is  the  shortwave  radiative  flux,  Qiw  is  the  longwave  radiative  flux,  Qs  is 
the  sensible  heat  flux,  and  Qi  is  the  latent  heat  flux.  All  fluxes  are  in  the  units  of 
watts  per  meter  squared  (W/m2).  These  terms  were  defined  as  positive 
downwards  and  vice  versa  as  archived,  and  thus  the  terms  were  simply  added 
together  to  yield  Qo. 

In  order  to  obtain  the  value  of  water  surface  friction  velocity,  uw*  the  u 
and  v-components  of  surface  wind  stress  were  first  combined  into  a  scalar  value 
of  surface  wind  stress,  in  units  of  N/m2.    This  value  of  wind  stress  was  then 
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converted  to  wind  speed  in  the  planetary  boundary  layer,  Vb,  in  units  of  meters 

per  second  (m/s).  The  following  equation  was  then  used  to  calculate  uw*: 

Uw*=(CDpa0/pw0)l/2\Vb\  (27) 

The  only  remaining  variable  in  the  calculation  of  P  is  h.  The  only 
modification  required  here  was  to  convert  the  depth  from  units  of  feet  to  meters. 
However,  since  the  grid  spacing  for  mixed  layer  depth  was  irregular  and  greater 
than  that  of  winds  and  heat  flux,  a  bivariable  interpolation  routine  was  used  to 
interpolate  the  polar  stereographic  mixed  layer  depth  grid  field  to  the  1.125 
Mercator  grid.  Contoured  plots  of  monthly  mean  P-values  are  shown  in 
Appendix  E. 

2.      Correlation  Computation 

The  first  step  in  calculating  the  correlation  between  dTjdt  and  h'  was  to 
match  the  grids  between  Ts  and  h.  For  ease  of  plotting,  the  2.5°  grid  for  Ts  was 
chosen.  Here  again  a  bivarible  interpolation  routine  was  used  to  interpolate  the 
polar  stereographic  mixed  layer  depth  grid  field  to  the  spherical  Ts  grid  field.  The 
objective  of  the  data  processing  here  is  to  obtain  a  correlation  field  based  on  daily 
trends  of  9T's/3t  and  h'  over  each  of  the  24  months.  Contoured  plots  of  these 
correlations  are  shown  in  Appendix  F. 
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III.    SURFACE  FIELDS 

The  daily  surface  data  sets  were  averaged  over  each  of  the  24  months  of  the 
two  year  period  to  yield  monthly  mean  averages. 

A.    HEAT  BUDGET 

In  order  to  see  the  time  variation  of  net  surface  heat  flux,  we  plot  the  zonal 
values  of  net  surface  heat  flux  versus  time  for  latitudes  5N,  0,  and  5S  (figure  7). 

1.  Western  Equatorial  Pacific 

The  heat  budget  in  the  western  equatorial  Pacific  is  characterized  by  a 
downward  net  surface  heat  flux  throughout  1986-87.  The  highest  downward 
heat  flux  occurs  in  the  southwestern  equatorial  Pacific  during  the  period  of 
September  -November  86  (225  W/m2).  In  the  northwestern  equatorial  Pacific,  the 
downward  heat  flux  values  oscillate  between  strong  during  the  spring  and  early 
fall  to  weak  during  the  summer  and  winter.  In  the  southwestern  equatorial 
Pacific,  downward  heat  flux  values  range  between  high  in  the  early  spring  and 
late  fall  to  weak  during  the  summer  months. 

2.  Mid  Equatorial  Pacific 

The  heat  budget  in  the  mid-equatorial  Pacific  is  characterized  by  a 
predominantly  upward  net  surface  heat  flux  over  the  course  of  the  two  year 
study.  The  maximum  upward  heat  flux  occurs  during  Aug  86  (255.9  W/m2)  in  the 
southern  portion  of  the  mid  equatorial  Pacific.  North  of  the  equator,  heat  fluxes 
are  strongest  during  the  winter,  diminishing  elsewhere.  South  of  the  equator, 
fluxes  are  again  predominantly  strongest  in  the  winter,  but  proceeding  further 
south,  fluxes  increase  strongly  in  the  summer. 
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3.      Eastern  Equatorial  Pacific 

The  heat  budget  in  the  eastern  equatorial  Pacific  has  a  typically 
downward  net  surface  heat  flux  over  a  period  of  two  years.  The  maximum 
downward  heat  flux  occurs  during  Feb  87  (199.2  W/m2)  in  the  north.  There  is, 
however,  no  distinct  seasonal  pattern  in  the  temporal  fluctuation  of  heat  flux  in 
this  region. 

B.    SURFACE  FRICTION  VELOCITY 

The  water  surface  friction  velocity,  uw*  as  described  by  equation  (5)  will  be 
referred  to  as  wind  speed  since  it  is  proportional  to  wind  speed  for  a  constant 
drag  coefficient.  The  time  evolution  of  zonal  winds  are  plotted  in  the  same 
manner  as  net  surface  heat  flux  and  are  shown  in  figure  8.  Meridional  winds  are 
not  plotted  due  to  their  low  relative  correlation  with  sea  surface  temperature 
change  as  shown  by  McPhaden  and  Hayes  (1991). 

1.      Western  Equatorial  Pacific 

a.  North  Equatorial  Latitudes 

The  western  equatorial  Pacific  is  dominated  by  easterly  winds  over 
most  of  the  two  year  period.  The  maximum  easterly  wind  speed  in  the  north  (.26 
cm/s)  occur  in  December  1986.  The  trend  in  the  north  is  typically  strong  easterlies 
during  the  winter  and  early  spring. 

b.  South  Equatorial  Latitudes 

The  maximum  easterly  winds  in  the  south  (.23  cm/s)  occur  in  June 
1986.  The  trend  in  the  south  is  typically  strong  easterlies  during  the  summer 
months  diminishing  to  weak  during  the  winter  months. 
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Figure  8.  Time  evolution  of  net  surface  heat  flux  (Q0)  for  latitudes  (a) 
5N,  (b)  0,  and  (c)  5S  across  the  equatorial  Pacific.  The  value  of  Q0  is 
defined  as  positive  (solid)  downwards,  and  negative  (dashed)  upwards. 
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2.  Mid  Equatorial  Pacific 

a.  North  Equatorial  Latitudes 

The  mid  equatorial  Pacific  is  a  region  of  seasonal  wind  variability. 
In  the  north  the  easterly  winds  dominate  in  late  winter  to  early  spring  with  the 
maximum  wind  speed  (.36  cm/s)  occurring  in  January  1987,  and  the  westerlies 
dominate  in  the  fall  with  maximum  wind  speed  (.38  cm/s)  occurring  in  September 
1986.  The  summer  months  represent  the  transition  between  the  easterlies  and 
westerlies  and  thus  the  near-zero  wind  speeds  are  fairly  consistent  with  the 
occurrence  of  the  doldrums. 

b.  South  Equatorial  Latitudes 

In  the  south,  the  seasonal  shifts  between  easterlies  and  westerlies, 
differ  from  those  in  the  north.  Easterlies  dominate  in  the  early  spring  to  late 
summer,  with  the  maximum  value  (.26  cm/s)  occurring  in  May  1987.  The 
transition  between  easterlies  and  westerlies  occur  here  in  the  early  fall  and  late 
spring.  The  westerlies  dominate  during  the  late  fall  and  winter  months,  with 
maximum  value  (.28  cm/s)  occurring  during  March  1987. 

3.  Eastern  Equatorial  Pacific 

The  eastern  equatorial  Pacific  follows  the  general  patterns  and  trends  of 
the  mid  equatorial  Pacific.  In  the  north,  however,  the  shift  between  easterlies  and 
westerlies  occurs  earlier,  in  the  spring,  and  thus  the  maximum  westerlies  (.61  cm/s) 
occur  in  the  summer,  specifically  July  1986.  The  maximum  easterlies  (.35  cm/s) 
occur  in  the  early  winter,  specifically  December  1986.  In  the  south,  the  winds  are 
fairly  weak  over  the  the  two  year  period. 
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Figure  9.  Time  evolution  of  water  surface  friction  velocity  (uw*)  for 
latitudes  (a)  5N,  (b)  0,  and  (c)  5S  across  the  equatorial  Pacific.  The 
value  of  uw*  is  defined  as  positive  (solid)  for  westerly  winds  and  negative 
(dashed)  for  easterly  winds. 
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C.    SEA  SURFACE  TEMPERATURE 

The  time  evolution  of  sea  surface  temperature,  Ts,  is  plotted  using  the  time- 
longitude  plots  used  in  the  previous  two  sections  and  is  shown  in  figure  9. 

1.  Western  Equatorial  Pacific 

The  western  equatorial  pacific  possesses  the  highest  values  of  sea 
surface  temperature  throughout  the  equatorial  Pacific.  These  high  values 
typically  increase  to  the  south.  The  highest  values  of  sea  surface  temperature 
occur  during  the  period  October- December  1986  and  March-July  87  (in  excess  of 
30Q  C)  in  the  southern  portion  of  the  mid-western  equatorial  Pacific.  Temporally, 
the  western  equatorial  Pacific  has  by  far  the  lowest  sea  surface  temperature 
variability  (  1Q  C).  These  findings  are  consistent  with  values  obtained  from 
Reynolds  NMC/CAC  Surface  Temperature  Analysis  and  the  findings  of 
McPhaden  and  Hayes  (1991)  and  McPhaden  (1988). 

2.  Eastern  Equatorial  Pacific 

The  eastern  equatorial  Pacific  has  the  lowest  values  of  sea  surface 
temperature,  generally  decreasing  to  the  southeast.  The  lowest  values  occur 
during  September  1986  and  August- September  1987  (about  20°  C).  The  general 
trend  of  sea  surface  temperature  shows  warm  values  propagating  eastward, 
peaking  in  March-April  of  each  year,  and  propagating  westward,  peaking  in 
September-October  of  each  year.  The  temporal  variability  of  temperature 
anomaly  was  also  plotted  versus  time  in  the  same  manner  as  the  latter.  The  results 
show  that  the  eastern  equatorial  Pacific  has  the  highest  variability  in  sea  surface 
temperature  [0(1  C)],  with  the  highest  variability  occurring  in  the  far  southeastern 
equatorial  Pacific  in  March  1986  (3°  C)  and  March  1987  (5Q  C). 
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Figure  10.  Time  evolution  of  sea  surface  temperature  (Ts)  for 
latitudes  (a)  5N,     (b)  0,  and  (c)  5S  across  the  equatorial  Pacific. 
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IV.    TWO  REGIMES  OF  THE  OML 

The  surface  forcing  function,  or  P-value  is  an  indication  of  the  dominant 
regime  of  the  ocean  mixed  layer,  whether  entrainment  (deepening),  or  shallowing. 
As  mentioned  previously,  the  surface  forcing  impacting  these  regimes  are  wind 
shear,  or  shear  production,  and  heat  flux,  or  buoyancy  damping.  These  fields 
have  been  described  in  the  chapter  III.  In  addition  to  these  variables,  the  mixed 
layer  depth  is  also  used  in  the  calculation  of  the  P-value. 

A.    MIXED  LAYER  DEPTH 

Mixed  layer  depth  is  plotted  and  analyzed  as  per  the  previous  sections.  The 
plot  of  time  evolution  of  mixed  layer  depth  versus  longitude  is  shown  in  figure 
10. 

L      Western  Equatorial  Pacific 

The  western  equatorial  Pacific  has  an  intermediate  range  of  mixed  layer 
values,  with  higher  values  increasing  to  the  east.  The  northern  region  has  the 
largest  values  of  mixed  layer  depth  during  the  spring  of  each  year  (30  m).  The 
lowest  values  of  mixed  layer  depth  occur  during  late  summer  each  year  (20  m) 
with  values  increasing  toward  the  equator.  The  southern  region  has  the  largest 
values  of  mixed  layer  depth  during  late  fall  of  each  year  (30  m),  with  the  lowest 
values  occurring  during  the  spring  of  each  year  (17-19  m). 

2.      Mid  Equatorial  Pacific 

The  western  portion  of  the  mid-equatorial  Pacific  has  largest  values  of 
mixed  layer  depth.  The  highest  values  of  mixed  layer  depth  occur  just  below  the 
equator  (50  m)  and  are  sustained  throughout  the  two  year  period  with  a 
comparatively  low  temporal  variability  (3-4  m). 
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3.      Eastern  Equatorial  Pacific 

The  eastern  equatorial  Pacific  maintains  the  lowest  values  of  mixed  layer 
depth,  particularly  in  the  far  eastern  portion.  The  lowest  values  (5  m)  occur  in  the 
north  (7.5  N)  during  the  spring  of  each  year.  The  lower  latitudes  (south  of  the 
equator)  tend  to  have  the  highest  values  of  mixed  layer  depth.  The  temporal 
trend  reveals  a  deepening  mixed  layer  in  the  fall  and  winter,  returning  to  a 
shallow  mixed  layer  depth  during  the  spring  and  summer.  This  trend  is  more 
pronounced  in  the  region  south  of  the  equator. 

B.    P-VALUE 

Based  on  equation  (5),  the  value  of  the  surface  forcing  function,  or  P-value,  is 
computed  using  the  water  surface  friction  velocity,  uw*,  the  mixed  layer  depth,  h, 
and  the  net  surface  heat  flux,  Qo.  A  plot  of  P-value,  averaged  over  two  years  at 
each  grid  point,  versus  time  is  shown  in  figure  1 1  for  latitude  5N,  the  equator,  and 
5S. 

1.  Western  Equatorial  Pacific 

The  western  equatorial  Pacific  is  dominated  by  negative  P-values.  In 
terms  of  the  mixed  layer,  this  would  tend  to  indicate  a  shoaling  regime.  The 
largest  negative  value  (2.3  (cm/s)3)  occurs  to  the  far  west  (150  E)  at  the  equator  in 
February  1986.  Although  there  is  some  temporal  variability  in  the  P-value  in  this 
region,  there  is  no  marked  trend. 

2.  Mid  Equatorial  Pacific 

The  mid  equatorial  Pacific  is  characterized  by  positive  P-values.  In  terms 
of  the  mixed  layer,  this  would  indicate  an  entrainment  regime.  The  largest 
positive  value  (2.49  (cm/s)3)  occurs  in  the  middle  (150  W)  just  south  of  the 
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Figure  11.  Time  evolution  of  mixed  layer  depth  (h)  at  latitudes  (a)  5N, 
(b)  0,  and  (c)  5N  across  the  equatorial  Pacific. 
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equator  in  August  1986.    There  is  no  temporal  trend  in  the  magnitude  of  the  P- 
value. 

3.      Eastern  Equatorial  Pacific 

The  eastern  equatorial  Pacific  has  both  positive  and  negative  P-values. 
The  positive  values  of  P  occur  toward  the  west,  blending  in  with  the  positive 
values  from  the  mid  equatorial  Pacific.  The  negative  values  occur  to  the 
southeast.  The  largest  negative  value  (-0.58  (cm/s)3)  occurs  to  the  far  south  in 
December  1986. 
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Figure  12.  Average  P-values  for  latitudes  (a)  ION  through  0,  and  (b)  0 
through  10S  across  the  equatorial  Pacific. 
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V.     CORRELATIONS  BETWEEN  3T's/dT  AND  H' 

Correlations  between  3T's/9t  and  h'  are  computed  with  h'  leading  3Ts/8t  by 
one  day.  Why  is  this  lag  time  used?  McPhaden  and  Hayes  (1991)  showed  that 
the  strongest  correlation  between  sea  surface  temperature  and  wind  forcing 
occurred  with  wind  forcing  leading  temperature  by  one  day  (figure  11).  It  can 
therefore  be  assumed  that  mixed  layer  depth,  which  is  driven  by  turbulence 
created  by  wind  mixing,  must  therefore  lead  temperature  on  the  same  time  scale. 
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Figure  13.  Cross  correlations  between  SST  and  wind  speed,  zonal  wind 
pseudostress,  and  wind  pseudowork.  Cross  correlation  extrema  (r)  and 
the  lag  (in  days)  at  which  they  occur  are  shown.  In  each  case,  the  lag 
indicates  that  cold  SST  follows  high  winds  by  1  day  (from  McPhaden  and 
Hayes,  1991). 
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In  order  to  study  this  correlation  and  its  relationship  to  P-value,  the  following 
process  was  used.  First,  dTJdt  and  h'  were  correlated  based  on  daily  values  over 
the  course  of  each  of  the  24  months,  yielding  monthly  correlation  fields.  These 
fields,  shown  in  Appendix  F,  are  shaded  dark  for  regions  having  negative 
correlations  and  light  for  regions  having  positive  correlations.  Comparing  these 
fields  to  the  P-flelds  shown  in  Appendix  E,  it  is  clear  on  the  whole  that  there  is  no 
similarity  between  respective  positive  and  negative  regions  as  would  be  expected 
by  "switching"  theory. 

The  emphasis  now  shifts  to  choosing  regions  in  the  equatorial  Pacific  that 
exhibit  typically  strong  entrainment  and  shallowing  tendencies,  based  upon 
positive  and  negative  P-values  respectively.  The  regions  chosen,  based  upon 
figure  11,  are  latitudes  (10N-10S)  along  165W  for  entrainment,  and  latitiudes 
along  160E  for  shallowing.  The  monthly  correlations  at  these  points  are  plotted 
over  time,  and  are  shown  in  figures  12  and  13.  The  result  here  is  that  over  time, 
both  of  these  regions  undergo  positive  and  negative  correlations  between  dTJdt 
and  h',  regardless  of  the  region  chosen  for  its  entrainment/shallowing  tendency.  It 
is  clear  from  the  apparently  random  oscillatory  behavior  of  these  correlations  that 
the  correlations  between  dTJdt  and  h'  that  we  would  expect  from  the  modified 
Kraus-Turner  model  are  not  dominant  on  time  scales  of  a  day  or  more. 

It  is  possible,  however,  that  although  the  correlations  between  dTJdt  and  h' 
are  not  significant  on  these  time  scales,  that  the  correlations  between  simply  Ts 
and  h'  might  be  significant.  A  test  of  this  hypothesis  is  meaningful  since  the 
correlation  between  these  two  variables,  T's  and  h'  is  also  commonly  accepted  as 
being  only  positive.  Correlations  between  daily  values  of  T's  and  h'  are  computed 
over  the  two  year  period  at  5S  and  longitudes  160E  and  90W  (shoaling  regions) 
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and  165W  (entrainment  region)  with  the  results  plotted  in  figure  16.  At  90W,  the 
highest  correlation  (-0.8)  occurs  with  h'  leading  T's  by  several  days.  The 
difference  at  this  longitude  is  that  instead  of  a  negative  maximum  correlation  for 
h'  leading  Ts,  the  maximum  correlation  is  positive  (.35)  and  occurs  with  h'  leading 
T's  by  nearly  100  days.  Of  importance  here  is  that  negative  correlations  are 
expected  between  temperature  and  mixed  layer  depth  in  regions  of  negative  P- 
values  (90W  and  165E)  and  conversely  positive  correlations  in  regions  of 
positive  P- values  (160E). 
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Figure  14.  Monthly  correlation  of  3T's/9t  and  h'  versus  time  for 
longitude  165W  and  (a)  latitudes  ION  -  0,  (b)  latitudes  0-10S. 
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Figure  15.  Monthly  correlation  or9T's/3t  and  h'  versus  time  for 
longitude  160E  and  (a)  latitudes  ION  -  0,  (b)  latitudes  0-10S. 
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Figure  16.  Correlation  between  daily  values  of  T's  and  h'  for  the  period 
1986-87  with  normalized  values  of  correlation  plotted  against  the 
lag(days)  at  latitude  5S  and  longitudes  (a)  90W,  (b)165 W,  and  (c)160E. 
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VI.    SUMMARY  AND  DISCUSSION 

In  studying  the  correlation  between  the  time  rate  of  change  of  sea  surface 
temperature  anomaly,  dTJdt,  and  the  mixed  layer  depth  anomaly,  h',  several 
findings  are  apparent.  The  most  significant  finding  is  that  there  exist  both 
positive  and  negative  correlations  between  dTJdt  and  h'  at  the  various  grid 
points  chosen  across  the  equatorial  Pacific.  This  finding  runs  contrary  to  the 
upwelling  physics  theory  showing  a  positive-only  correlation.  Since  forward 
differencing  was  used  in  calculating  dTJdt,  the  correlations  were  based  on  a  one 
day  lead  of  h'.  The  reason  for  using  this  lead  time  is  that  mixed  layer  depth  and 
sea  surface  temperature  vary  with  mesoscale  atmospheric  variability  (Elsberry  and 
Garwood,  1978),  and  mixed  layer  depth  changes  caused  by  wind  stress  induced 
turbulence  lead  sea  surface  temperature  changes  by  about  one  day  (McPhaden 
and  Hayes,  1991). 

A  second  finding  is  that  in  the  western  equatorial  Pacific  (140E-180E)  the 
monthly  averaged  values  of  the  surface  forcing  function,  P,  are  typically  negative 
throughout  the  course  of  the  two  year  period,  indicating  that  this  region  is 
predominantly  a  shoaling  regime.  This  is  due  largely  to  the  strong  downward  net 
surface  heat  flux  in  this  region.  In  the  mid-equatorial  Pacific  (180E-120W),  the 
monthly  averaged  values  of  the  surface  forcing  function,  P,  are  typically  positive 
indicating  that  this  region  is  predominantly  an  entrainment  regime.  Finally,  the 
eastern  equatorial  Pacific  is  characterized  by  both  positive  and  negative  values, 
with  negative  values  typically  to  the  east  indicating  a  general  shallowing  regime. 

Although  the  correlation  patterns  between  3T's/3t  and  h'  do  not  fit  the  P- 
value  distribution  with  any  degree  of  confidence,  there  is  some  similarity  between 
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the  correlations  of  Ts  and  h'  and  P-values.  It  is  found  that  h'  leads  Ts  by  several 
days  with  a  strong  negative  correlations  (-0.8,  -0.6)  in  regions  of  significantly 
negative  P-values  (90W,  160E).  The  converse  is  not  so  conclusive  in  that  h'  leads 
T's  by  nearly  100  days  with  a  relatively  strong  positive  correlation  (0.35)  in  a 
region  of  positive  P-values  (165W). 
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VII.     CONCLUSIONS  AND  RECOMMENDATIONS 

A  positive-only  correlation  between  dTJdt  and  h'  is  commonly  accepted  in 
current  thermodynamic  theories  regarding  the  El  Nino  phenomenon.  In  order  for 
these  theories  to  be  valid,  this  relationship  must  hold  true,  particularly  in  the 
equatorial  Pacific.  Based  on  the  results  of  this  study,  this  positive-only 
correlation  is  not  true  in  the  equatorial  Pacific.  It  would  be  premature  at  this 
point,  though,  to  discount  the  validity  of  these  theories  based  solely  on  the  results 
of  this  study.  To  more  sufficiently  study  this  correlation  and  its  relation  to  the  P- 
value,  the  following  study  is  recommended.  First,  the  equatorial  Pacific  should  be 
broken  up  into  three  regions,  western,  mid,  and  eastern,  to  investigate  the  unique 
surface  forcings  occurring  in  each.  In  each  of  these  regions,  thermistor  chains  or 
current  meters  should  be  utilized  at,  above  (about  5N),  and  below  (about  5S)  the 
equator  to  measure  hourly  surface  and  subsurface  temperatures.  The  subsurface 
temperatures  could  then  be  analyzed  to  determine  mixed  layer  depth  anomaly 
and  subsequently  correlated  with  the  time  rate  of  change  of  sea  surface 
temperature  anomaly.  This  finer  time  scale  would  not  only  conform  more  closely 
to  the  time  scale  of  turbulence  in  the  mixed  layer,  but  would  allow  more  flexibility 
in  adjusting  the  lead  time  used  in  the  correlations  between  mixed  layer  depth  and 
sea  surface  temperature.  The  conflicting  theories  regarding  this  correlation  could 
then  be  better  resolved. 
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